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ABSTRACT  . 
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F-irsri^measureivients  of  the  time  for  the  recovery  of  the  critical 

7^1  tv/ 

current  at  a phase  slip  site  have  been  made  bv  studyino  tho 
frequency  dependence  of  the  hysterotic  rf  current-voltage  charac- 
teristics in  the  low  power  limit.  The  data  are  explicable  by  a 
simple  quasiparticle  diffusion  model  and  are  incompatible  with 
recent  hot  spot  models. 


The  production  and  decay  of  nou'-cqu i 1 ibr ium  states  in 
superconductors  is  currently  of  great  interest.  In  recent 
studies,  non-equilibrium  states  have  been  produced  by 
tunneling  , laser  excitation  , thermal  phonon  injection  , 
ultrasonic  phonons^,  conversion  of  normal  to  super currents 
at  superconducting  boundaries  , rapid  variation  of  the 
supercurrent^ , etc.  In  this  paper  we  report  the  first 
determination  of  the  temporal  response  of  the  non-equilibri- 
um state  generated  by  the  phase  slip  process  at  a site  in  an 
ultra-thin-f ilm  superconducting  microbridge  in  a rf  SQUID 
configuration.  The  results-  are  inconsistent  with  the  hot 

7 

spot  model  but  agree  v/ell  with  a modified  quasiparticle 
8 

diffusion  model  . 

Our  data  were  obtained  using  a broadband  continuounljr 
variable  nonresonant  mutual  inductance  bridge  {lOMHz  - 
1.8GHz),  a conventional  20MHz  tuned  SQUID  biasing  circuit, 
and  a 9.2  GHz  spectrometer  employing  a resonant  cavity. 

In  order  to  minimize  heating  effects,  we  investigated 
numerous  thin  film  cylindrical  SQUIDs  at  20  MHz  and  9.2  GHz 
and  selected  a very  low  critical  current,  nearly  granular 
niobium  sample.  This  sample  allowed  us  to  study  the  hyster- 
etic  response  over  a wide  temperature  range  with  minimal 
heating.  The  physical  dimensions  of  this  niobium  wea)«  link 
are  as  follows:  3 micrometers  in  the  direction  of  current 
flow,  by  40  micrometers  wide,  by  35  Angstroms  thick,  totally 
covered  by  a 600  Angstrom  insulating  oxide  layer.  This 
sample  has  a superconducting  critical  temperature  of  2.8  K. 


No  chaiKjcs  in  tho  rf  rcr.ponse  wore  observed  with  the  sample 
craersed  in  liquid  helium  just  al.)ovc  or  belov;  tho  helium  x 
point  or  in  a vacuum  can. 

Insert  (a)  of  Fig.  1 illustrates  the  detected  rf  I—V 
characteristics  of  the  sample  for  various  frequencies  at 
2.167  K as  observed  with  the  broadband  system.  The  tempera- 
ture dependence  of  the  critical  current  can  be  described  irr 
terras  of  a mean  field  behavior.  This  behavior  is  indicative 
of  a current  distribution  v/hich  is  relatively  uniform  over 
the  width  of  the  raicrobridge.  Above  a frequency  of  120  UTlz 
hysteresis  develops  in  the  critical  currents.  That  is^ 
dissipation  persists  in  the  weak  link  belov/  the  meart  field 
critical  current  until  a sudden  return  occurs  (dov;n  arrov.'s) 
to  the  uniform  superconducting  state.  Note  that  the  mean 
field  critical  current  (up  arrows)  is  independent  of  fre- 
quency. - . 

This  critical  current  hysteresis  is  similar  to  that 

g 

observed  by  SBT  in  the  dc  I-V  characteristics  of  tin 
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microbr idges  and  by  Rachford  et  al.  in  the  9.2  GHz 
microwave  response  of  weak  links  of  various  geometries  and 
materials.  This  type  of  behavior  follows^®  from  the 
assumption  that  and  are  given  by: 

. ItI  V(C  (T)/3)  , (la) 

. ItI^  v/2a,  (lb) 

where  is  the  mean  field  critical  current  of  the  micro- 
bridge,  Y is  the  superconducting  ordorparameter , 5 (T)  is  the 
coherence  length,  and  is  the  minimum  current  necessary 
to  support  Josephson  oscillations  across  a distance  corapac- 


able  to  the  inelastic  scattering  length  The  length  a is 

1/2 

given  by  the  eguation,  A - (iv„t/3)  ' where  t is  the  elastic 

r 

length,  is  the  Fermi  velocity,  and  j is  the  appropri- 
ate inelastic  scattering  time.  The  ratio  of  the  squares  of 
the  critical  currents  is  proportional  to  temperature; 

= [2A/(n(0) /3)  ]^(T-T^)/T^.  (2) 

To  explain  the  frequency  dependence  of  the  hysteresis 
[illustrated  in  Fig.  1,  insert  (a)],  consider  an  increasing 

V 

instantaneous  rf  current  through  the  v;eak  link.  As  the 
magnitude  of  the  rf  current  exceeds  the  mean  field  critical 
current,  a phase  slip  site  is  rapidly  established  producing 
non-equilibrium  quasiparticles  and,  hence,  dissipation  in 
the  sample.  The  quasiparticles  depress  the  critical  current 
so  that  at  some  point  in  the  decreasing  current  part  of  the 
rf  cycle  (determined  by  a)  the  current  becomes  less  than 
this  new  cricital  value.  When  the  magnitude  of  the  current 
is  less  than  the  depressed  critical  current,  the  voltage 
across  the  site  vanishes,  the  excess  quasiparticle  popula- 
tion is  no  longer  sustained  by  the  current  and  begins  to 
decay.  The  critical  current  recovers  towards  its  mean  field 
value.  On  the  next  half  cycle  this  time  dependent  critical 
current  is  exceeded,  phase  slip  occurs,  and  the  non-equi- 
librium state  is  regenerated.  If  the  rf  period  is  long 
compared  to  the  decay  time  [see  Fig.  2,  insert  (a) 1 , the 
critical  current  will  recover  to  the  mean  field  value  and  no 
time  average  lowering  of  the  critical  current  (no  hystere- 
sis) will  be  observed. 


Hov.’over,  it  the  rf  period  is  r.liort  coir.parocj  to  the 
decay  time,  the  excess  qvias  ipar  1 5 clc-s  will,  not  dociiy 
Gufficionily  oetwe-en  halt  cycles.  In  this  Cci.e,  the  ciitical 
current  v/ill  not  liave  time  to  recover  to  the  ;:can  £i.cld 
value,  and  there  will  be  a net  lov.’ering  of  the  detected 
critical  current  as  long  as  the  non-cqu  ilibr irm  state  is  . 
cyclically  sustained.  This  condition  v.’ill  be  maintained 
until  the  rf  current  is  reduced  below  I^,(C)  [see  Fig.  2 
insert  (b)  ] . Thus  the  hysteresis  deperids  on  the  relative 
values  of  tlie  rf  period  and  the  decay  time  for  the  excess 
quasiparticle  distribution. 

i 

In  figure  1 , the  measured  values  of  [I  /r.,(f)l 

V t 

are  plotted  as  a function  of  temperature  at  various  frequen- 
cies. Note:  (1)  that  at  each  frequency  the  current  ratio 
can  be  represented  by  straight  linos  extrapolating  to  a 
single  (critical)  temperature;  (2)  there  is  a temperature 
dependent  frequency  for  the  onset  of  the  hysteresis;  C3) 

that  the  linear  slopes  approach  a constant  lii.iting  value  at 
higher  frequencies.  According  to  Eq.  2,  the  linear  beliavior 
implies  that  all  the  temperature  dependence  is  contained  in 
the  coherence  length  and  that  a is  temperaturr  i ndependent 
over  the  range  studied.  The  variation  of  the  irlopos  follows 
from  the  frequency  dependence  of  I^(f)  described  above. 


In  order  to  account  for  the  frequency  dcpondonco  of 

I ^(f) , it  is  sufficient  to  make  a simple  extension  of  the 

quasiparticlo  diffusion  model.  Assume  that  the  depressed 

critical  current  recovers  as  exp(  + t/T^P  or,  equivalently, 

that  the  quasiparticle  diffusion  distance  a decays  in  time 

as  A = Aq  exp(-t/T^),  where  is  a temperature  independent 

time  constant.  Also  assume  that  the  phase  slip  center  is 

established  in  a time  short  in  comparison  to  the  decay  time. 

The  minimum  rf  amplitude,  U (f)],  necessary  to  cyclically 

sustain  the  dissipative  state  at  the  phase  slip  site  occurs 

when  the  critical  current  (grovjing  in  time  as  a relaxes) 

tangentially  meets  the  instantaneous  rf  current  [see  Fig.  2, 

insert  (b) . ] Using  this  condition  and  the  above  assumptions, 

t^  [sec  Fig. 2,  inserts  (a)  and  (b) ] can  be  iteratively 

computed  using  ^^2  “ ~ i^^ln  (sin  wtj^/sin  (jt2)  / where  tj  is 

given  by  ~ ( tan~^  ut^)  / u.  Having  found  t2,  the  value  of 

I^(f)/I^  can  easily  be  obtained,  where  I ^ is  row  taken  to 

represent  the  saturated,  high  frequency  value  of  the  reduced 

o 

critical  current,  Usinq  this  rest.  . ind  (lyi.)  '»  the  ratio 
j[l|/I  Jr}]  is  obtained. 

The  solid  curve  in  Fig.  2 represents  a fit  of  this 
model  to  the  experimental  frequency  variation  of  (I  ^I^(f)I 
(solid  circles)  as  determined  by  the  slopes  of  the  lines 
shown  in  Fig, 2.  Two  parameters  we/e  used  in  obtaining  the 
solid  curve:  the  ratio  [I  and  the  constant 

The  shape  of  the  curve  is  determined  from  the  model  itself 
and  the  assumption  of  exponential  decay  of  the  non-oqui- 
libtium  state.  We  find  a saturated  value  a 3.87, 


9.0. 


and  a tempetaturo  independent  time  constant  = 5.5 
nanoseconds,  and,  therefore,  by  equation  (2),  ^‘q/Cq  = 

Estimating  ^(0)  from  the  bulk  by  the  formula 
?(0)  -■  find  = 2590  /vngstroins,  and  from 

A = (*-v,,t/3)  we  find  a quasiparticle  inelastic  scatter- 

r 

ing  time  t = 8.5  xlO  sec.  We  expect  that  x sets  the  time 
scale  for  the  establishment  of  the  quasiparticle  distribu- 
tion, \7hereas  is  the  time  for  the  relaxation  of  the 
non-equilibrium  state  involving  the  interaction  of  quasi- 
particles, 2A  phonons,  and  pairs^^. 

From  the  dissipation  observed  in  the  saraplc  at  9 GHz  we 

can  estimate  the  resistive  length  in  the  saniple  independent 

of  the  hysteresis.  Previously  it  has  been  observed  that  the 

resistive  length  is  approximately  twice  the  quasiparticle 

diffusion  length  A.  In  insert  (b)  of  Fig.  1 v;e  plot  the 

R 

quasiparticle  diffusion  length,  A , obtained  from  dissipa- 
tion  measurements.  Note  that  A is  independent  of  temper- 
ature and  is  in  excellent  agreement  v/ith  the  value  of  A 
obtained  from  the  development  of  hysteresis. 

Recently,  hysteresis  in  the  dc  I-V  characteristics 

of  weak  links  has  been  explained  in  terras  of  a simple  hot 
7 12 

spot  model  ' . A straightforward  application  of  this 

model  predicts  that  the  temperature  difference  AT  between 
the  bath  and  the  hotspot  should  decay  as  1 - e~^^^ . 

A more  detailed  analysis  of  the  heat  balance  equation. 


v;hero 


there  in  heating  on.Ty  if  the  instantaneous 


tenperature 


exceeds  T or  if  the  rf  current  exceeds  the  temperature 
c 

depcndeiit  critical  current,  still  results  in  an  asymptotic 
relaxation  of  &T-  Any  asymptotic  relaxation  of  at  and, 
hence,  the  local  critical  current  to  the  mean  field  value, 
will  require  that  the  temperature  at  which  hysteresis  is 
first  observed  is  independent  of  frequency.  This  result 
is  inconsistent  v^ith  our  observations  (sec  Fig.l.). 

In  the  low  pov;er  regime,  we  mfiy  expect  that  quasipar- 
tides  v:ill  be  trapped  by  Andreev-like  scattering  iu  the 
vicinity  of  a phase  slip  site  due  to  a local  reduction  of 
the  average  superconducting  gap  in  this  region.  The  scale 
of  the  region  of  trapping  will  be  set  by  inelastic  scatter- 
ing processes  of  the  quasiparticles  and  will  sJirink  to  zero 
as  the  excess  local  excitations  relax.  If  this  is  true, 
then  A can  be  identified  v.’ith  the  width  of  the  quasi- 
particlo  trapping  region  v;hich  relaxes  in  a characteristic 
time  T^. 

In  our  experiments,  the  phase  slip  process  is  synchron- 
ized with  the  rf  cycle?  whereas,  in  dc  I-V  experiments, 
phase  slip  proceeds  at  the  Josephson  frequency.  In  most 
cases  the  Josephson  period  is  fast  v;ith  respect  and  the 
observed  hysteresis  is  consistent  with  the  usual  analysis 
with  a temperature  independent  A.  However,  if  the  critical 
current  is  small  and  if  the  voltage  across  the  site  produces 
Josephson  oscillations  with  a period  of  the  order  of 


then  a voltuige  dcpenclcncc  of  A wiJl  be  r,oc?n.  TiiiiJ  effect 
simply  icflocts  the  f requency  (t ime)  clcpcndencc  of  A cliscuGscd 
above. 

.In  summary,  we  have  presonted  the  first  i.iea5.:uremont 
of  the  broadband  temporal  response  of  a "weal:  .link".  Our 
results  have  been  interpreted  over  the  entire  temperature 
range  by  a simple  extension  of  the  qunsiparticlc  diffusion 
model.  The  data  rire  inconsistent  v;iLh  simple  liotspot  theory 
or  any  model  which  predicts  an  asymptotic  recovery  of  the 
critical  current  to  the  mean  field  value.  Froi.i  our  measure- 
ments we  have  determined  times  characteristic  -of  the 
dy  narnics  of  the  non-cquilibr itm  state  produced  by  the  phase 
slip  process. 
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Fig.  1.  ' ohtainocl  from  data  such  ac  shov/n  in 

insert  (a)  ic  plotted  as  a function  of  temperature.  Insert 
(a)  shows  the  observed  development  of  iiystercsis  in  the 
detected  I-V  curves  as  frequency  is  increased.  Insert  (b) 
shows  the  temperature  independence  of  the  half-length,  , 
of  the  dissipative  region  as  infered  from  the  GlI/: 
dissipation  mcasureiuonts. 

Fig.  2.  The  slopes,  d [I^/I^  (f ) ] VdT,  from  tlie  previous  figure 
are  plotted  as  a function  of  frequency  (solid  circles).  The 
solid  curve  is  calculated  from  the  model  described  in  the 
text  v.’ith  = 5.5  nsec,  and  (I^/I^)^  = 3.87.  In  inserts  (a.) 

and  (b)  , t^^  is  the  time  at  which  the  magnitude  of  the  rf 
current  (solid  line)  exceeds  the  time  dependent  critical 
current  (dashed  line)  and  t2  is  the  time  at  v.’hich  the  rf 
current  is  less  than  I.rtn<(.i  (f)  is  the  lov;eGt  (frequency 

v V 

dependent)  current  that  will  cyclically  sustain  the  dissipative 


state  shown  in  (b) . 
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